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REAL-TIME TOOL FOR AFFINE TRANSFORMATIONS OF 
TWO DIMENSIONAL IFS FRACTALS 


ELENA HADZIEVA AND MARIJA SHUMINOSKA 

Abstract. This work introduces a novel tool for interactive, real-time 
transformations of two dimensional IFS fractals. We assign barycentric 
coordinates (relative to an arbitrary affine basis of R 2 ) to the points that 
constitute the image of a fractal. The tool uses some of the nice prop¬ 
erties of the barycentric coordinates, enabling any affine transformation 
of the basis, done by click-and-drag, to be immediately followed by the 
same affine transformation of the IFS fractal attractor. In order to have 
a better control over the fractal, as affine basis we use a kind of minimal 
simplex that contains the attractor. We give theoretical grounds of the 
tool and then the software application. 


1. Introduction 

One of the simplest ways of creating fractals is by means of iterated func¬ 
tion systems (IFSs). Generally, the IFS can be defined on any complete 
metric space, but for the purposes of this paper we will restrict our discus¬ 
sion to the two-dimensional real space R 2 . In this section we will go through 
the basic definitions, given in mm- The transformation defined on R 2 of 
the form 

w(x) = A • x + b, 

where A is 2 x 2 real matrix and b is a two-dimensional real vector, is called 
a (two-dimensional) affine transformation. The finite set of affine contrac¬ 
tive transformations, {r;i,r; 2 , ... ,w n } with respective contractivity factors 
i = 1, 2,..., n, together with the Euclidean space (R 2 , d^;) is called (hy- 
perbolyc) iterated function system (IFS). Its notation is {R 2 ; wi,W2, ..., w n } 
and its contractivity factor is s = max{s^,z = 1,2, Since (R 2 ,^) 

is a complete metric space, such is (%(R 2 ),/i(d^;)), where (^(R 2 ), d(d#)) is 
the space of nonempty compact subsets of R 2 , with the Hausdorff metric 
h(dE ) derived from the Euclidean metric dE (0)- It can be shown that 
if the hyperbolic IFS {R 2 ; wi, 1^2, • • •, w n} with the contractivity factor s is 
given, then the transformation W defined on 7d(R 2 ), by 

W(B) = U i = 1 Wi(B), B G 
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is a contraction mapping with the contractivity factor s. According to the 
fixed-point theorem, there exist a unique fixed point of W, A G H(M 2 ), 
called the attractor of the IFS, that obeys W(A) — A — lim^^oo W(B ), for 
any B G (Note that there are weaker conditions under which the 

attractor of an IFS can be obtained - the IFS need not to be contractive, 
see |3].) 

Two algorithms are used for the visualization of an IFS attractor, deter¬ 
ministic and random (the last is also called chaos game algorithm). In this 
paper we construct the images of the attractors with the second, more effi¬ 
cient algorithm. Our primary target in this paper will be fractal attractors, 
since there is no interactive, real-time tool for their modeling, to the best of 
our knowledge.. 

The rest of the paper is organized in four sections: Related Articles, 
Theoretical Grounds of the Tool, The Tool and its Application, Conclusions 
and Future Work. 


2. Related Articles 

Barnsley et ah, in 0 E] define the fractal transformation as mapping 
from one attractor to another, using the “top” addresses of the points of 
both attractors and illustrate the application in digital imaging. Although 
not user friendly, not real-time and not continuous (they are continuous 
under certain conditions), these transformations, have a lot of potential for 
diverse applications. In [7] the coordinates of the points that compose the 
fractal image are determined from the IFS code and then the IFS code is 
modified to obtain translated, rotated, scaled or sheared fractal attractor, 
or attractor transformed by any affine transformation as the composition of 
aforementioned transformations. To make a new transformation, another 
IFS code has to be constructed. In [8], Darmanto et al. show a method 
for making weaving effects on tree-like fractals. They make a local control 
of the “branches” of the tree, by changing a part of the IFS code. The 
control over the change of the attractor depends on the experience of the 
user for predicting how the modification of the coefficients in the IFS code 
will change the tree-like fractal. Compared to the methods for transforming 
fractals in the cited papers, our tool is user-friendly, real time, relatively 
fast and relatively low memory consuming. It enables a continuous affine 
transformation of an arbitrary IFS attractor. 

The idea for modeling fractals, i.e. their predictive, continuous transfor¬ 
mation, where barycentric coordinates are involved, is exposed in mm, 
[Tj, by means of, so-called Affine Invariant Iterated Function Systems - 
AIFS. The IFS code is transformed into AIFS code which involves barycen¬ 
tric coordinates. The method that we propose is also based on barycentric 
coordinates, but no transformation of the IFS code is needed. 

Also, Kocic et al. in M propose the use of, so called minimal canonical 
simplex , for better control of the attractor (for 2D case, minimal canonical 
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simplex is the isosceles right-angled triangle with the minimal area that 
contains the attractor and whose catheti are parallel to the x and y axes). 
In )2\ the authors prove a theorem for existence and uniqueness of such 
minimal canonical simplex. It is the limit of the Cauchy sequence of the 
minimal simplexes of the n-th preattractor when n approaches to infinity, in 
the complete metric space ( / H(R m ), /i(d^)). (The n-th preattractor in this 
paper is defined as the set obtained after n iterations of the chaos game 
algorithm.) 


3. Theoretical Grounds of the Tool 

Barycentric coordinates of points (and vectors) promote global geomet¬ 
rical form and position rather than exact coordinates in relation to the 
origin/axes. Thus, the geometric design done by the use of barycentric 
coordinates can be called coordinate-free geometric design. 

The set of the three points { A , £>, C} is said to be affine basis of the affine 
space of points R 2 if the set {A^, A ( 3 } is a vector base of R 2 , observed as 
a vector space (mi)- We say that the point M has barycentric coordinates 
(a, 6, c) relative to the basis {A , £>, C}, where a + b + c = 1, and we write 

M = aA + bB + cC (1) 

if and only if one of the following three equivalent conditions holds: 

Til = bA^ + cAp 

M = obX + cB(5 
cPi = acl + bci 


We will give the relations between the barycentric coordinates relative to 
the given basis {A,B,C} and the rectangular coordinates of an arbitrary 
point from R 2 . Suppose that the arbitrary point M and the three basis 
points A, B and C have rectangular coordinates (x,?/), (^ 1 ,^ 2 )? (& 1 ,^ 2 ) an d 
(ci, C 2 ), respectively. Then the relation 0 can be rewritten in the following 
form 


or, in matrix form 
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By rearranging the last matrix equation, we obtain more suitable form 
for further manipulations: 
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The relation ^ defines the conversion of the barycentric coordinates rel¬ 
ative to the affine basis {A, B,C} into the rectangular coordinates. Let us 
denote the 3x3 matrix by T to indicate “triangle”. The inverse conversion 
can be easily get, by simple multiplication of <© by the inverse of T (which 
exists, since {A, B,C} is an affine basis). That is, the conversion of the 
rectangular coordinates into the barycentric coordinates relative to the basis 
{A , B , C} is defined by the relation 


a 


X 

b 

c 

1—1 

I 

H 

II 

1 

1 _ 


( 4 ) 


which after calculating the inverse of T, can be expressed in its explicit form, 


a 

b 
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i>2 — c 2 ci - bi bic 2 - b 2 c\ 


X 

C2 — CL2 a\ — ci a 2 ci — aic 2 


y 

_ a 2 -b 2 b\- a\ a\b 2 - a 2 &i 


1 


( 5 ) 


where det T = a±b2 — b\a2 + b\C2 — C 1&2 + 01*22 — aiC 2 ^ 0. 

Barycentric coordinates have two other names: affine and relative coor¬ 
dinates. Affine, since they are related to an affine basis, and relative, since 
they define relative positions in the plane/space. Namely, if we consider a set 
of points with given barycentric coordinates relative to the basis {A, £>, C} 
and we affinly transform the basis into the basis {A , ,S / ,C / }, then the set 
of points having the same barycentric coordinates, but now relative to the 
basis {A 7 , 5', C"}, will keep the relative geometry, i.e. will be transformed 
by the same affine transformation. Note that any change of one affine basis 
to another, defines a unique affine transformation. 

Since the image of the attractor, generated by the random algorithm 
is a finite set of points, we will use the aforementioned property of the 
barycentric coordinates to affinly transform the IFS fractal. 


4. The Tool and Its Application 

Our real time, user-friendly tool is written in Cff programming language 
using the Visual Studio 2013. It allows both on-click definition of the triangle 
and definition by specifying the triangle’s vertices directly in the code, with 
their rectangular coordinates. When the affine base is defined, the image of 
the fractal is created (relative to the triangle) and ready for transformation. 

Triangle’s vertices are moved by drag and drop option of the cursor. When 
the cursor is placed over the vertices, it has a “hand” shape, otherwise its 
shape is “cross”. The affine transformations of the triangle, immediately 
followed by the same affine transformation of the fractal image, are visible 
in real time, in the coordinate system of the window. 

We will show the superiority of our tool on two examples of fractal at¬ 
tractors, the so called “flower” and “maple”. In order to get clearer image 
of the attractor, we neglect the first 14 points and start to plot from the 
15-th point. 
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Note that, by default, the origin of the coordinate system of the window 
is located at the top-left corner of the window, with positive directions: to 
the right for the x-axis and down for the y-axis. 

Example 1. The IFS {M 2 ; wi, whose attractor is the “flower” fractal, 
is defined by the following contractive mappings : x i-A A^x + b^, i = 1,2: 

0.37 " 

1.74 J 5 

-0.34 " 

1.75 * 

Figure [I] (a), b), c)) depicts that the arbitrary control triangle is defined 
on-click. When moving its vertexes by click-and-drag, the fractal attractor 
appropriately responds to the transformations done over the triangle. 

Example 2. The attractor of the IFS {M 2 ; w i, IT 2 , ^3, ^ 4 }, where the map¬ 
pings : x i-A A^x + b i, i = 1, 2, 3,4, are defined by 

" -0.08 " 

0.26 J 5 

0.07 " 

3.5 J 5 

0.74 " 

0.39 J 5 

-0.56 " 

0.60 J 5 

We computed the coordinates of the vertices of the minimal canonical 
simplex for this fractal and used the vertices as an affine basis. Such affine 
basis ensures better control over the attractor an, m , see Figure [2] We 
examined creation of the shadow or wind effect for the maple tree, through¬ 
out d) to f) on this figure. 

5. Conclusions and future work 

We succeeded in creating a real-time and user-friendly tool for interactive 
modeling of IFS attractors, with focus on fractal attractors. To the best 
of our knowledge, this is the first real-time interactive tool for modeling 
fractals. The tool is very efficient because it is relatively low memory and 
time consuming. We foresee a great application of our tool, since all images 
coded by an IFS code can be affinly transformed according the user needs, 
only by click-and-drag. 

These are the first steps in fractal modeling with our tool. There is much 
more to be done to enhance the tool’s performances. For instance, including 
a control polygon of more than three points. Moreover, it would be very 
efficient to use the convex hull of the attractor as a control polygon. In the 
last case the control over the global changes of the attractor will be maximal. 
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Including local control would be also useful for satisfying the sophisticated 
needs of the users. Finally, expanding the tool over 3D IFS attractors will 
bring higher practical importance to the tool. 

References 

[1] E. Babace: Iterated Function System with Affine Invariance Property (in Macedo¬ 
nian), PhD Thesis, Faculty of Natural Sciences and Mathematics, University ”Ss. 
Cyril and Methodius”, Skopje, R. Macedonia, 2009 

[2] E. Babace, Lj. M. Kocic : Minimal Simplex for IFS Fractal Sets, NAA 2008, Lecture 
Notes in Computer Science 5434, pp.168-175, Eds.: S. Margenov, L. G. Vulkov, J. 
Wasniewski, Springer-Verlag Berlin Heidelberg, 2009 

[3] M. F. Barnsley: Fractals everywhere, Academic Press, 1988 

[4] M. F. Barnsley, A Vince: The Chaos Game on a General Iterated Function System, 
Ergodic Theory Dynam. Systems 31, no. 4, 1073-1079, 2011 

[5] M. F. Barnsley, B. Harding: Fractal Transformations in 2 and 3 Dimensions, 
arXiv:1308.6648vl, 30 Aug 2013 

[6] M. F. Barnsley, T. Harding, K. Igudesman, How to transform and Filter Images 
Using Iterated Functions Systems, SIAM Journal of Imaging Sciences,vol. 4, no. 4, 

pp 1001-1028, 2011 

[7] H. T. Chang: Arbitrary affine transformation and their composition effects for two- 
dimensional fractal sets, Image Vision Comput. 22 (13) , 1117-1127, 2004 

[8] T. Darmanto, I. S. Suwardi, R. Munir: Weaving Effects in Metamorphic Animation 
of Tree- like Fractal based on a Family of Multi-transitional Iterated Function System 
Code, Computer, DOI: 10.1109/IC3INA.2013.6819150 Conference: The International 
Conference on Computer, Control, Informatics and its Applications 2013 (IC3INA 
2013), At Jakarta, Indonesia, Volume: I 

[9] T. D. DeRose: Three-Dimensional Computer Graphics, A coordinate-Free Approach 
(manuscript), University of Washington, 1992 

[10] K. J. Falconer: Fractal Geometry. Mathematical foundations and Applications, John 
Wiley and Sons, 1990 

[11] J. Gallier: Geometric Methods and Applications For Computer Science and Engi¬ 
neering, Springer - Verlag, TAM Vol.38, 2000 

[12] Lj. M. Kocic: AIFS - a Tool for Biomorphic Fractal Modeling, Nonlinear Dynam., 
Psychology Life Sci. 5, no. 1, 45-63, 2000 

[13] Lj. M. Kocic, A. C. Simoncelli: Towards free-form fractal modelling, In: Mathemat¬ 
ical Methods for Curves and Surfaces II, M. Daehlen, T. Lyche and L. L. Schumaker 
(eds.), pp. 287-294, Vanderbilt University Press, Nashville (TN.), 1998 

[14] Lj. M. Kocic, L. Stefanovska, E. Babace: AIFS and the Minimal Simplex Problem, 
Proceedings of The International Conference of Differential Geometry and Dynamical 
Systems, 5-7.10.2007, Bucharest, Romania, pp.119-128, 2008 

University ”St. Paul the Apostle”, Partizanska bb, Ohrid, Macedonia 
E-mail address : elena.hadzieva@uist. edu.mk 


University ”St. Paul the Apostle”, Partizanska bb, Ohrid, Macedonia 
E-mail address : marija. shuminoska@cse .uist. edu.mk 



REAL-TIME TOOL FOR AFFINE TRANSFORMATIONS OF TWO DIMENSIONAL IFS FRACTALS 


o 

A + 



a) 


b) 



Figure 1. a) Drawing point A; b) Drawing point B; c) 
Drawing point C and the fractal; d) Moving point B, after 
point C is moved; e) Moving point A; f) Moving point C. 
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Figure 2. a) Creating simultaneously the minimal canonical 
simplex and the fractal relative to the simplex; b) Moving 
point B; c) Moving point C; d) Moving point A; e) Moving 
point C; f) Moving point C further. 













